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SECTION  I 
INTRODUCTION 

In  a  1984  study  the  Energy  Conversion  Branch  of  the  U.S.  Air  Force 
Wright  Aeronautical  Laboratories  Aero  Propulsion  Laboratory  investigated  the 
use  of  advanced  power  generating  devices  for  future  mobile  electric  power 
(MEP)  applications.  One  conclusion  of  that  study  was  that  kinematic 
Stirling  engines  have  a  high  potential  for  meeting  the  requirements  for  mid¬ 
sized  (30-100  kW)  flightline  and  operational  system  electronics  support 
applications. 

There  are  several  kinematic  Stirling  engines  currently  under  develop¬ 
ment  by  various  organizations.  An  example  is  the  automotive  Stirling  engine 
now  being  developed  under  the  Department  of  Energy  (DOE)  Automotive  Heat 
Engine  Program  with  program  management  by  NASA  Lewis  Research  Center.  One 
current  effort  also  sponsored  by  this  DOE  program  is  Contract  DEN3-351  with 
Stirling  Thermal  Motors,  Inc.  (STM)  of  Ann  Arbor,  Michigan,  to  experi¬ 
mentally  evaluate  advanced  Stirling  concepts  [1].  STM,  using  primarily 
private  funding,  has  designed,  and  is  fabricating,  a  Base  Technology 
Stirling  Engine,  designated  the  STM4-120  (4  cylinders,  120  cc  swept  volume 
per  cylinder),  incorporating  these  advanced  concepts.  The  STM4-120  design 
is  suitable  for  a  variety  of  applications  and  since  its  heat  input  is  from  a 
liquid  metal  heat  pipe  it  can  be  connected  easily  to  almost  any  heat  source. 

The  STM4-120  design  features  several  advanced  concepts  which  should 


reduce  the  Stirling  engine  size,  weight,  complexity  and  expected  manufac¬ 
turing  costs.  Indirect  heating  technology,  incorporating  heat  pipes,  is  an 
integral  part  of  the  STM4-120.  making  it  possible  to  divide  the  engine  into 


an  energy  conversion  unit  (ECU)  and  a  distinctly  separate  external  heating 
system  (EHS),  thus  permitting  simplification  of  the  heat  exchanger  design. 

The  heater  of  the  ECU,  actually  the  uniform  temperature  condenser  of  the  EHS 
heat  pipe,  is  designed  to  take  advantage  of  the  high  film  coefficient  of  the 
condensing  metal  vapor.  This  permits  design  optimization  based  on  Stirling 
cycle  thermodynamic  requirements,  and  not  on  flue  gas  heat  transfer, 
resulting  in  improved  engine  performance. 

Other  advanced  features,  such  as  the  variable  angle  swashplate  power 
control  and  compliantly  mounted  reciprocating  seals,  are  presented  in  the 
literature  [2]. 

An  analysis  by  Argonne  National  Laboratory  looked  specifically  at  the 
STM4-120's  Air  Force  MEP  applicability  and  concluded  that  such  use  is 
attractive  [3]. 

These  two  studies  and  the  STM4-120  design,  led  the  Air  Force  Wright 
Aeronautical  Laboratories  to  sponsor  an  additional  task  to  the  current  DOE 
sponsored  activity  at  STM.  The  new  Air  Force  sponsored  task  was  to  design 
the  external  heating  system  for  the  STM4-120,  which  would  meet  the  appro¬ 
priate  military  specifications,  MIL-STD-633  and  MIL-G-52884  for  a  30  kW 
mobile  electric  power  generator  set. 

The  STM4-120  was  designed  to  produce  about  40  kW  (53  HP)  at  3000  RPM. 
The  Air  Force  has  an  MEP  requirement  for  a  range  of  20  to  60  kW.  Although 
the  STM4-120  will  produce  approximately  25  kW  (33  HP)  at  1800  RPM,  the 
heating  system  will  be  designed  to  the  full  rated  power  of  40  kW.  The 
eventual  reoptimization  of  the  STM4-120  to  produce  the  military  required 
power  at  1800  RPM  is  a  straightforward  procedure  and  will  have  relatively 
little  impact  on  the  external  heating  system  design.  A  complete  description 


of  the  STM4-120  is  included  in  Appendix  A. 

The  Air  Force  Wright  Aeronautical  Laboratories  requested  NASA  Lewis 
Research  Center  to  undertake  program  management  responsibility  through 
Military  Interdepartmental  Purchase  Request  Number  FY  1455-84-N0618,  to 
expand  the  effort  at  STM  to  include  the  conceptual  design  of  an  external 
heating  system  and  the  preliminary  design  of  a  heat  pipe  heater  for  a 
Stirling  powered  electrical  generator  set. 
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SECTION  II 

EXTERNAL  HEATING  SYSTEM  CONCEPTUAL  DESIGN 

1.  EXTERNAL  HEATING  SYSTEM  CONFIGURATION 

The  conceptual  design  of  the  external  heating  system  (EHS)  was  based  on 
the  heat  input  requirements  of  the  STM4-120  Base  Technology  Stirling  Engine. 
The  EHS  was  designed  to  operate  on  liquid  fuel,  using  a  flat  plate,  counter 
flow  recuperator  and  a  corrugated,  flat-surface  heat  exchanger.  In  both 
components  the  flow  of  air  and  flue  gas  is  strictly  laminar  (Poiseuille  flow 
between  flat  plates)  and  hence  the  Nusselt  number  is  independent  of  the 
Reynolds  number.  The  Nusselt  number  based  on  gap  hydraulic  diameter,  is 
taken  to  be  always  7.6  [4]. 

The  combustion  chamber  design  was  chosen  from  an  existing  Philips 
Stirling  engine,  the  ADVENCO  4-88.  The  burner,  atomizing  air  nozzle  and 
mixing  chamber  were  readily  adaptable  to  the  EHS  design  concept  regarding 
fuel,  size  and  performance. 

In  operation,  preheated  air  is  mixed  with  atomized  fuel  in  the  combus¬ 
tion  chamber  and  burns  continuously.  The  atomized  fuel  is  introduced  into 
the  combustion  chamber  through  a  nozzle.  The  geometry  of  the  combustion 
chamber  was  designed  for  proper  mixing  and  distribution  of  flue  gas  flow 
into  the  heat  exchanger.  Primary  air  is  fed  into  the  combustion  chamber 
through  a  swirl  chamber  and  secondary  air  is  fed  downstream  through  radial 
holes.  Film  cooling  further  downstream  prevents  overheating  of  the 
combustion  chamber  walls. 

Hot  flue  gas  flows  from  the  combustion  chamber  through  the  multiple 
parallel  gaps  of  the  heat  exchanger,  transferring  its  heat  through  the  wall 
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into  the  finned  enclosure.  The  enclosure  is  the  evaporator  section  of  a 
heat  pipe  where  liquid  sodium  in  the  wick  lining  the  inside  wall  evaporates. 

The  sodium  vapor  flows  to  the  engine  via  connecting  tubes. 

Since  the  evaporator  enclosure  was  designed  to  operate  at  sub-atmospheric 
pressure,  coarse  wire  mesh  was  incorporated  to  provide  structural  support  to 
the  fin  walls,  while  providing  a  vapor  path  for  the  sodium  working  fluid. 

The  residual  heat  of  the  flue  gas  leaving  the  heat  exchanger  is  used  to 
preheat  the  incoming  combustion  air  in  the  recuperator. 

A  schematic  diagram  of  the  EHS  is  shown  in  Figure  1.  Due  to  the  use  of 
heat  pipes  for  heat  transport  the  EHS  can  be  located  a  variety  of  locations 
relative  to  the  STM4-120  ECU. 

An  analysis  and  simulation  code  for  the  EHS  was  developed  to  determine 
the  influence  of  various  geometric  parameters  on  the  system  performance. 

The  range  of  parameter  variation  was  limited  by  the  following 
considerations: 

•  Packaging  requirements  limited  the  overall  system  dimensions  to 
about  330  mm  x  250  mm  x  100  mm  as  shown  in  Figure  2. 

•  Heat  exchanger  (evaporator)  fins  and  preheater  plates  were 
designed  with  the  same  depth  for  ease  of  fabrication  and  packaging 
convenience. 

•  In  laminar  flow  heat  exchangers  the  Nusselt  number  can  be  assumed 
to  be  independent  of  the  Reynolds  number,  therefore  it  was 
advantageous  to  divide  the  flow  path  into  as  many  parallel  paths 

as  possible,  to  derive  the  maximum  reduction  in  flow  friction 
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Figure  1  -  External  Heating  System  Schematic  Illustration 


(1)  Fuel  atomizing  nozzle  (2)  Combustion  chamber  (3)  Heat  exchanger/evaporator 
(4)  Air  preheater  (5)  Exhaust  (6)  Supporting  wire  mesh  (7)  Fine  gauze  wick 
(8)  Evaporator  fin 
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power  loss,  as  shown  in  Table  1.  The  narrowest  gap  deemed 
practical  was  1.0  mm. 

2.  SYSTEM  ANALYSIS  AND  OPTIMIZATION 

The  geometrical  parameters  of  the  recuperator  (or  preheater)  and  heat 
exchanger  alone  determine  the  performance  of  the  EHS  provided  that  proper 
mixing  and  residence  time  can  be  achieved  in  the  combustion  chamber.  Since 
an  existing  design  of  proven  satisfactory  performance  was  specified  for  the 
burner,  the  combustion  chamber  design  was  decoupled  from  the  system  analysis 
and  optimization. 

a.  Description  of  the  System  Analysis  Code 

The  purpose  of  the  analysis  and  simulation  code  was  to  determine 
the  EHS  efficiency  and  the  effect  on  flow  losses  and  cold  start  fuel 
consumption  based  on  the  geometric  parameters  of  the  recuperator  and  heat 
exchanger. 

Aerodynamic  and  heat  transfer  models  were  employed  to  define  the  design 
parameters  based  upon  the  heat  transfer  between  the  exhaust  and  intake  air 
in  the  recuperator,  the  combustion  characteristics  of  the  fuel  and  burner, 
and  the  heat  transfer  from  the  flue  gas  to  the  heat  pipe  evaporator  [5]  [6]. 

In  heat  transfer  calculations,  the  temperature  dependent  properties  of 
the  air  and  flue  gas  were  taken  into  consideration  [7]  and  therefore,  a 
finite  differences  numerical  scheme  was  used  to  solve  the  energy  equations. 
Furthermore,  in  calculation  of  the  flame  temperature  and  heat  exchanger 
performance,  the  dissociation  of  the  flue  gas  was  taken  into  consideration. 


This  effect  was  neglected  in  the  analysis  of  the  recuperator  since  the 
temperatures  are  sufficiently  low. 

Since  the  total  laminar  pressure  drop  in  the  heat  exchangers  was 
expected  to  be  less  than  50  cm  of  water  it  was  deemed  unnecessary  to  take 
into  account  the  variable  fluid  properties  in  the  flow  direction.  The 
pressure  drop  was  determined  for  each  component  separately,  based  on 
Poiseuille  flow  at  the  average  temperature  prevalent  in  that  component. 

These  average  temperatures  were  calculated  by  numerical  integration. 

The  code  was  designed  to  accommodate  any  fuel  (atomic  hydrogen  to 
carbon  ratios  between  0.25  and  4.0)  and  equivalence  ratios  ranging  from  1.0 
to  1.4. 

The  governing  equations  for  the  preheater  and  heat  exchanger/evaporator 
heat  transfer  analyses  are  included  in  Appendix  B.  A  sample  computer 
run  is  included  in  Appendix  C. 

b.  Results  of  the  System  Analysis  and  Optimization 

In  order  to  investigate  the  influence  of  the  ratio  of  the  evapora¬ 
tor  area  to  the  recuperator  area  (with  a  constant  sum  thereof)  a  series  of 
simulations  was  made  with  the  following  geometric  parameters  unchanged; 
fin/plate  depth,  number  of  gaps  and  gap  width  in  both  the  evaporator  and 
recuperator,  and  the  sum  of  the  lengths  of  the  evaporator  and  recuperator. 

The  results  are  summarized  in  Figures  3  and  4  showing,  respectively, 
the  effect  of  the  area  ratio  on  the  efficiency  and  the  cold  start  fuel 
consumption.  The  area  ratio  is  defined  as  the  ratio  of  recuperator  area  to 
heat  exchanger/evaporator  area.  The  cold  start  fuel  consumption  decreases 
with  an  increase  of  the  area  ratio,  and  efficiency  exhibits  a  maximum  at  an 


area  ratio  value  of  approximately  0.25.  A  design  point  of  0.32  was  chosen 
to  account  for  possible  deviation  of  the  analytical  model  from  reality  and 
to  allow  sufficient  interface  space.  That  area  ratio  was  held  constant  for 
the  rest  of  the  investigation. 

The  effect  of  total  heat  transfer  area  at  a  constant  ratio  of  heat 
transfer  area  to  cross  sectional  area  was  investigated  by  a  series  of 
simulations  with  increasing  depth  of  the  evaporator  fins  and  recuperator 
plates  maintaining  all  other  geometric  parameters  constant.  The  results  are 
summarized  in  Figures  5,  6,  and  7,  showing  respectively,  the  effect  of  the 
total  heat  transfer  area  on  the  efficiency,  flow  losses  and  cold  start  fuel 
consumption.  A  design  point  at  the  total  area  of  4.235  m  corresponding  to 
a  fin/plate  depth  of  70  mm  was  selected  as  a  favorable  compromise  between 
efficiency  and  flow  losses  on  one  hand  and  cold  start  fuel  consumption 
(size,  weight,  and  cost)  on  the  other  hand. 

It  should  be  noted  that  the  efficiency  and  cold  start  fuel  consumption 
in  Figures  3  and  4  differ  from  the  values  in  Figures  5  and  7.  This  is  due 
to  the  iterative  procedure  used  to  determine  first  the  design  point  with 
regard  to  area  ratio,  then  the  design  point  with  regard  to  total  heat 
transfer  area.  The  final  value  of  total  area  is  larger  than  the  value  used 
to  determine  area  ratio. 

Finally,  the  sensitivity  of  the  performance  to  the  gap  width  was 
investigated  by  a  series  of  simulations  with  varying  gap  width.  The  number 
of  gaps  for  each  simulation  was  determined  based  on  the  gap  width  so  that 
the  total  width  of  the  heat  exchangers  was  maintained  constant  at  250  mm. 
The  results  of  this  investigation,  summarized  in  Figures  8  and  9,  shows  as 
expected,  that  increasing  the  gap  width  has  an  adverse  effect  on  the 


efficiency  which  is  much  more  severe  than  the  beneficial  effect  on  flow 
losses.  Consequently  a  gap  width  of  1.0  mm  was  retained  as  the  design 
point. 

The  entire  investigation  was  performed  at  full  load  corresponding  to 
DF2  fuel  consumption  of  2.67  grams/sec  with  an  equivalence  ratio  of  1.25 
(air/fuel  ratio  of  18:1),  4%  atomizing  air  and  no  EGR.  The  material  wall 
thickness  was  kept  at  0.2  mm.  The  goals  of  this  investigation  were  to 
establish  preliminary  design  parameters,  hence  the  effects  of  fouling  and 
corrosion  were  not  taken  into  consideration. 

The  performance  at  the  above  conditions  with  the  selected  geometric 
parameters  is  summarized  in  Table  2  and  Figures  10a  and  10b. 
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Table  2 


STM4-120  EHS  Performance  Data 


Preheater 


Evaporator 


Geometric  Data 

Air  Side 

Flue  Side 

Length,  mm 

220 

220 

110 

Depth,  mm 

70 

70 

70 

Gap  width,  mm 

1.0 

1.0 

1.0 

Number  of  gaps 

104 

104 

67 

Wall  thickness,  mm 

0.2 

0.2 

0.2 

Temperature  data 

Average  (bulk)  temperatures,  K 

694 

S27 

1370 

Fuel  Data  (DF2) 

Atomic  H/C  ratio 

1.804 

Specific  lower  heating  value 

10.179  cal./gr. 

Operating  Data 

Equivalence  ratio 

1.25 

Atomizing  air 

4% 

Fuel  consumption 

2,67  gr./sec. 

Ambient  temperature 

301  K 

Heat  pipe  temperature 

1085  K 

Performance 

Preheated  air  temperature 

1034  K 

Flame  temperature 

2397  K 

Heat  exchanger  outlet  temperature 

1097  K 

Exhaust  outlet  temperature 

498  K 

Heat  output 

104  kW 

EHS  Efficiency 

91.6% 

Cold  start  penalty 

128.3  gr. 

Pressure  drop 

4.64  k  Pa 

Flow  losses 

819  Watt 
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SECTION  in 


HEAT  PIPE  PRELIMINARY  DESIGN 

The  conceptual  Stirling  powered  electrical  generator  set  includes  the 
STM4-120  Energy  Conversion  Unit  (ECU)  coupled  to  the  EHS  through  four  heat 
pipe  tubes  with  mechanical  couplings  that  simultaneously  supply  the  ECU  with 
heat  via  sodium  vapor  and  return  the  sodium  condensate  to  the  EHS. 

In  this  section,  a  brief  review  of  current  heat  pipe  theory  is 
presented,  followed  by  a  description  of  the  EHS  heat  pipe  configuration  and 
the  application  of  the  theory  to  the  particular  design  for  the  EHS. 

1.  Heat  Pipe  Theory 

The  theory  of  heat  pipes  is  well  documented  in  the  literature,  hence 
only  a  brief  review  of  current  theory  is  presented  to  provide  a  basis  for 
the  hydrodynamic  and  thermal  design  of  the  EHS  heat  pipe.  Parts  of  this 
section  are  excerpted  from  "Theory  and  Design  of  Variable  Conductance  Heat 
Pipes,"  by  D.B.  Marcus,  NASA  Report  Number  CR-2018.  A  listing  of  other 
pertinent  references  is  included  in  that  report. 

a.  Hydrodynamics 

Because  a  heat  pipe  involves  the  circulation  of  a  working  fluid, 
certain  pressure  drops  arise.  In  general,  there  will  be  viscous  losses  due 
to  liquid  flow  in  the  wick  or  capillary  structure,  and  viscous  and  inertial 
losses  due  to  vapor  flow  in  the  core.  In  addition,  there  may  be  body  forces 
that  either  aid  or  hinder  circulation  (e.g.  acceleration  fields  due  to 


etc.).  For  steady  state  operation  of  a  heat  pipe,  a  pressure  head  equal  to  the 
sum  of  these  losses  must  be  supplied  by  the  capillary  in  the  wick.  This  yields 
the  following  steady  state  pressure  balance,  which  must  be  satisfied  between 
all  points  along  the  heat  pipe. 


AP  =  AP„  +  AP  ±  AP. 
c  a  v  b 


net 

capillary 

head 


liquid 

pressure 

drop 


+ 


vapor 

pressure 

drop 


body  force 
head 
(if  any) 


(1)  Capillary  Head 

Of  the  terms  in  Eq.  1,  the  liquid  and  vapor  pressure  drops  are 
functions  of  the  circulation  rate  and  increase  with  the  heat  transfer  load  on 
the  device,  while  the  body  force  term  is  usually  independent  of  the  load.  Thus, 
to  satisfy  Eq.  1,  the  capillary  head  must  also  increase  with  load  in  such  a  way 
as  to  match  the  losses  incurred. 

The  capillary  head  in  a  saturated  heat  pipe  wick  arises  as  a  dynamic 
phenomenon.  It  is  due  to  the  existence  of  a  pressure  difference  across  a  curved 
liquid-vapor  interface  which  is  given  by: 


APi  =  o  (- 


where  APi  =  Interfacial  pressure  difference 


o  =  Surface  Tension 


Rj,  r  Two  orthogonal  radii  of  curvature  of  the  interface 


This  pressure  difference,  which  for  concave  menisci  results  in  a  depres¬ 
sion  of  the  liquid  pressure  with  respect  to  the  vapor,  exists  all  along  the 
heat  pipe  wick.  In  order  to  obtain  a  net  capillary  head,  it  is  necessary  for  it 
to  be  greater  at  the  evaporator  than  at  the  condenser.  In  a  heat  pipe  under 
load,  this  is  exactly  what  occurs  due  to  changes  in  the  interface  curvatures. 
Vaporization  of  the  liquid  in  the  evaporator  causes  the  menisci  to  recede  into 
the  wick  resulting  in  a  decrease  in  the  radii  of  curvature,  while  condensation 
in  the  condenser  has  the  opposite  effect.  Therefore,  capillary  pressure  is  not 
constant.  Capillary  pumping  is  a  passive  phenomenon  which  automatically  adjusts 
to  meet  the  flow  requirements,  within  limits. 

The  extremes  of  the  interfacial  pressure  difference,  which  can  be  gener¬ 
ated  along  the  wick  depend  on  the  wick  geometry  and  the  wetting  angle  between 
the  liquid  and  wick  material.  For  the  EHS  application,  the  wick  was  specified 
to  be  stainless  steel  wire  mesh,  based  on  previous  experience.  For  this  type  of 
wick,  the  maximum  capillary  pressure  becomes: 

P  =  (3) 

c  r 

max  c 

where  rc  =  effective  capillary  radius,  =  w  *  — 
w  =  wire  spacing 
d  =  wire  diameter 

(2)  Liquid  Pressure  Drop 


Liquid  flow  in  heat  pipes  is  generally  characterized  by  very 


low  velocities  and  Reynolds  numbers.  Consequently,  inertial  effects  can  be 
neglected  for  steady  state  operation  and  the  flow  losses  attributed  only  to 
viscous  shear. 

The  pressure  drop  in  wire  mesh  wicks  can  be  expressed  as  a  form  of  Darcy's 
Law  for  flow  in  porous  media,  which  is  usually  expressed  in  terms  of  a  permea¬ 
bility  -  K,  a  measure  of  the  wick  flow  resistance. 


In  Eq.  5,  d  is  the  wire  diameter  and  e  can  be  calculated  by  the 


equation: 


e  =  1  -  ir/4  SNd  (6) 

where  N  =  mesh  number 

S  =  crimping  factor  si. 05 

(3)  Vapor  Pressure  Drop 

The  vapor  pressure  drop  in  heat  pipes  is  often  considerably 
more  difficult  to  calculate  than  that  in  the  liquid,  for  in  addition  to 
viscous  shear,  the  analysis  must  account  for  momentum  effects  and  perhaps 
turbulent  flow  and  compressibility.  Complicating  this  analysis  is  the  fact 
that  mass  addition  in  the  evaporator  and  mass  removal  in  the  condenser  can 
significantly  alter  the  velocity  profiles  and  hence,  the  local  pressure 
gradient. 

From  conservation  of  momentum,  the  vapor  pressure  drop  can  be  expressed 
as  follows: 


dP 
_ v 

dx 


-(f  Re  )  u  m 

V  V  rV  V 

2  A  r  2  p 

V  V  V 


-e 


(7) 


Here,  the  first  group  represents  the  pressure  drop  due  to  viscous  shear 
and  the  second  group  represents  the  pressure  drop  due  to  inertial  effects, 
where  f  is  the  frictional  drag  coefficient  and  £  is  defined  by  the 
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equation: 


6  = 


/  V2  dA 


A 


v 


(8) 


Equation  8  accounts  for  the  effects  of  changes  in  vapor  velocity  across 
the  section  of  consideration. 

(4)  Body  Force  Head 

The  last  term  in  the  pressure  balance  equation  is  the 
pressure  head  due  to  body  forces  acting  on  the  working  fluid.  In  the  EHS 
heat  pipe,  the  only  body  force  that  arises  is  due  to  the  acceleration  of 
gravity.  Since  the  density  of  the  vapor  is  very  much  lower  than  that  of  the 
liquid,  the  body  forces  on  the  vapor  can  be  neglected.  The  pressure  head 
due  to  the  gravitational  field  is  expressed  by  the  equation: 


AP  =  J  dPK  =  ±J  p£  g  cosQdz  (9) 

°  dz 

where  AP  =  component  of  the  hydrostatic  head  in  the 
6 

direction  of  liquid  flow 
g  =  acceleration  of  gravity 

0  =  angle  of  the  liquid  flow  axis  with  respect  to 
the  direction  of  gravitational  acceleration 


dz  =  elemental  length 


.v.v.vyyLV.vv' 


(5)  Capillary  Pumping  Limit 

The  fact  that  there  exists  a  maximum  capillary  head  for  any 
wick-fluid  combination  results  in  a  hydrodynamic  limit  on  heat  pipe 
capacity.  As  mentioned  previously,  the  capillary  head  must  increase  with 
the  liquid  and  vapor  pressure  drops  as  the  heat  load  increases.  Since  there 
exists  a  limit  on  the  capillary  head  there  must  also  exist  a  corresponding 
limit  on  the  heat  load  if  the  pressure  balance  criterion  is  to  be  satisfied. 
This  defines  the  capillary  pumping  limit,  expressed  by  the  following 
inequality: 


AP  -  AP  S  AP^  +  APv  UO) 

cmax  ® 

(6)  Entrainment  Limit 

A  phenomenon  which  can  effect  the  capillary  pumping  limit, 
but  was  not  included  in  Eq.  10,  is  liquid  entrainment  in  the  vapor.  In  a 
heat  pipe,  the  vapor  and  liquid  generally  flow  in  different,  often  opposite, 
directions.  Since  they  are  in  contact  at  the  wick  surface,  this  sets  up  a 
mutual  drag  at  the  interface.  If  the  relative  velocity  between  the  liquid 
and  the  vapor  becomes  too  great,  the  interface  becomes  unstable  and  droplets 
of  liquid  will  be  torn  from  the  wick  and  entrained  in  the  vapor.  Since  this 
liquid  never  reaches  the  evaporator,  it  cannot  contribute  to  the  heat 
transferred  by  the  heat  pipe.  However,  it  does  contribute  to  the  liquid 
flow  loss.  Thus,  the  maximum  axial  heat  transfer  in  the  heat  pipe  is  no 
longer  equal  to  the  maximum  fluid  circulation  rate  times  the  latent  heat  of 
vaporization,  but  some  lower  value  which  defines  the  entrainment  limit. 
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The  conditions  leading  to  entrainment  are  expressed  in  terms  of  the 
ratio  of  vapor  inertial  forces  to  liquid  surface  tension  forces,  called  the 
Weber  number: 


p  V  2z 
v 


where  =  vapor  density 

V  =  average  vapor  velocity 


o  =  liquid  surface  tension 

z  =  characteristic  dimension  associated  with  liquid 


surface 


Limited  experimental  data  with  screen  wicks  indicate  that  a  Weber 
number  of  unity  represents  the  entrainment  condition  when  the  charac¬ 
teristic  dimension,  z,  is  set  approximately  to  the  screen  wire 
diameter  [8]. 

When  the  Weber  number  is  set  equal  to  unity,  the  limiting  axial 
heat  flux  corresponding  to  incipient  entrainment  is  given  by: 


p  a 

[  ~~z 


where  Q  =  Maximum  heat  flux  for  incipient  entrainment 
emax 

A  =  Vapor  core  flow  area 
v 

X  =  Latent  heat  of  vaporization 
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(7)  Sonic  Limitation 

It  can  be  shown  analytically  [9]  that  there  is  a  corres¬ 
pondence  between  constant  area  flow  in  a  heat  pipe  with  mass  addition 
(evaporator)  and  removal  (condenser)  and  constant  mass  flow  in  a  converging- 
diverging  nozzle.  The  end  of  the  evaporator  corresponds  to  the  throat  of 
the  nozzle.  Consequently,  just  as  there  is  a  sonic  limitation  on  the  flow 
velocity  through  a  nozzle  throat,  there  is  a  similar  limit  on  the  flow 
velocity  at  the  heat  pipe  evaporator  exit.  For  a  given  exit  temperature  and 
working  fluid,  this  choked  flow  condition  is  a  fundamental  limitation  on  the 
axial  heat  transfer  flux  capacity  of  the  heat  pipe. 

The  sonic  limit  is  calculated  by  setting  the  vapor  flow  velocity  equal 
to  the  sonic  velocity  in  the  continuity  equation  and  multiplying  by  the 
latent  heat  of  vaporization  as  follows: 


max 


=  X  p  V 
v  s 


(13) 


where  Os  =  axial  heat  flux  at  Mach  1  conditions 

'  max 

V  =  sonic  velocity  of  the  vapor 


b.  Heat  Transfer 

The  fluid  circulation  phenomenon  discussed  in  the  previous 
sections  arise  as  a  result  of  heat  transfer  into  the  heat  pipe  at  the 
evaporator  and  out  of  it  at  the  condenser.  In  most  heat  pipes  this  heat 
must  be  transferred  through  the  walls  and  saturated  wick.  Generally,  these 
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processes  are  the  major  source  of  temperature  drop  in  the  heat  pipe. 
Temperature  drops  also  arise  due  to  vapor  flow  losses  along  the  heat  pipe 
and  to  non-equilibrium  at  the  liquid-vapor  interfaces. 

(1)  Boiling  Limitation 

The  fact  that  heat  is  transferred  into  the  heat  pipe  through 
the  wick  gives  rise  to  another  limit  on  heat  pipe  capacity.  Frequently, 
liquid  is  vaporized  only  at  the  wick  surface  as  a  result  of  heat  conducted 
through  the  wick.  However,  the  vapor  at  the  wick  surface  is  thermo¬ 
dynamically  saturated,  hence  the  fluid  within  the  wick  is  superheated  by 
virtue  of  the  curvature  of  the  menisci  and  the  temperature  gradient  in  the 
saturated  wick.  The  greatest  superheat  occurs  at  the  interface  of  the  wick 
and  wall.  If  this  superheat  becomes  too  large  (it  increases  with  the  heat 
transfer  rate),  the  fluid  will  begin  to  boil  within  the  wick.  This  results 
in  evaporator  dryout  due  to  a  reduction  in  liquid  flow  area  in  the  wick 
caused  by  nucleation. 

Further  discussion  of  boiling  limitation  is  referred  to  the  literature, 
since  the  heat  flux  was  specified  to  be  a  value  experimentally  proven  to  be 
within  the  boiling  limitation  of  the  EHS  evaporator  configuration. 

(2)  Heat  Pipe  Temperature  Characteristics 

Generally,  heat  pipe  temperature  characteristics  are 
evaluated  at  full  load,  since  the  associated  gradients  are  at  a  maximum 
under  these  conditions.  The  temperature  drops  associated  with  heat  pipes 
occur  in  the  evaporator  wall,  the  saturated  wick,  and  the  condenser  wall  due 
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to  conduction.  There  is  an  additional,  generally  much  smaller,  temperature 
drop  associated  with  the  vapor  pressure  drop  between  the  evaporator  and  the 
condenser  due  to  vapor  flow  resistance. 


(a)  Temperature  Drop  Across  the  Evaporator  Wall 


The  temperature  of  the  outer  wall  of  the  evaporator  is 
essentially  constant  at  1085  K  (812°C).  The  temperature  drop  due  to  conduction 
across  the  wall  can  be  expressed  by  the  following  equation: 


_SL  - 

A  " 


(14) 


(b)  Temperature  Drop  Across  the  Saturated  Wick 

The  temperature  drop  across  the  saturated  wick  can  be 
evaluated  using  Eq.  14.  However,  the  thermal  conductivity  must  be  replaced 
by  an  effective  conductivity  of  the  wick/liquid  combination.  The  effective 
conductivity,  k  ,  for  saturated  screen  wicks  is: 
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where  k  ^  =  Liquid  Thermal  Conductivity 
kw  =  Wick  Thermal  Conductivity 
e  =  Wick  Porosity 
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(c)  Temperature  Drop  Due  to  Vapor  Flow  Resistance 


The  temperature  drop  of  the  vapor  due  to  flow  resistance 
can  be  evaluated  using  the  Clausius-Clapeyron  equation: 
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(16) 


(d)  Temperature  Drop  Across  the  Condenser  Tubes 

The  condensers  for  the  EHS  heat  pipe  are  the  STM4-120 
engine  heaters.  These  condensers  are  simple  tube  bundles  sized  to  optimize 
the  heat  transfer  to,  and  flow  resistance  of,  the  engine  working  fluid.  The 
temperature  drop  across  the  condenser  tubes  can  be  evaluated  by  modifying 
Eq.  14  to  cylindrical  coordinates.  The  resulting  equation  is: 
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2.  Heat  Exchanger/Evaporator  Configuration 

The  STM4-120  Energy  Conversion  Unit  (ECU)  is  coupled  to  the  EHS  through 
four  heat  pipe  tubes  with  mechanical  couplings  that  simultaneously  supply 
the  ECU  with  heat  via  sodium  vapor  and  return  the  sodium  condensate  to  the 
EHS.  The  heat  pipe  configuration  is  shown  schematically  in  Figure  11. 

Liquid  return  is  accomplished  via  gravity  by  locating  the  EHS  heat  exchanger 
below  the  ECU. 

Tests  conducted  at  STM  on  the  mechanical  couplings,  Aeroquip  Conoseals, 


Tube  containing  Stirling 
Engine  working  fluid 


Figure  11 

Schematic  Representation  of  STM4-120  Heat  Pipe 


available  from  the  Marman  Division  of  Aeroquip,  proved  the  functionality  and 
durability  of  these  seals  in  sodium  heat  pipe  applications.  One  heat  pipe 
with  this  coupling  has  operated  for  more  than  1775  thermal  cycles  of  heating 
to  800°C  and  cooling  to  150°C.  Another  heat  pipe  with  its  wick  capillary 
connection  across  this  coupling  has  operated  against  gravity  for  more  than 
1000  hours.  These  tests  are  continuing  [1].  An  exploded  view  of  a  Conoseal 
joint  is  shown  in  Figure  12. 

The  EHS  heat  pipe  configuration  was  examined  as  both  buffered  (moving 
front)  and  non-buffered  designs.  The  buffered  configuration  permits  easy 
start-up  and  temperature  control  but  suffers  from  the  fact  that  the  sharp 
temperature  gradients  in  the  front  create  extremely  high  thermal  stresses, 
especially  in  the  engine  heater.  As  the  front  passed  through  this  heat 
exchanger,  some  of  the  tubes  would  be  at  800°C  and  others  would  still  be  at 
ambient,  resulting  in  unacceptably  high  stresses,  and  increasing  the 
probability  of  heater  tube  failure. 

In  a  non-buffered  configuration,  the  entire  heat  pipe  increases  in 
temperature  uniformly,  eliminating  the  difficulties  that  result  from  thermal 
stresses.  However,  this  configuration  can  experience  start-up  difficulties 
by  condensing  and  freezing  the  sodium,  drying  out  the  evaporator  before  full 
fluid  circulation  can  be  established. 

The  EHS  heat  pipe  fluid  circulation  was  designed  to  return  the  maximum 
amount  of  sodium  to  the  evaporator  upon  shut-down  by  wicking  only  the 
evaporator.  This  technique  is  only  applicable  if  liquid  return  from  the 
condenser  can  be  accomplished  without  capillary  pumping.  In  the  EHS  heat 
pipe,  liquid  return  is  accomplished  via  gravity. 

Start-up  is  more  readily  accomplished  since  raising  the  temperature  of 


the  working  fluid  occurs  at  the  heat  source,  not  at  remote  sections  of  the 
heat  pipe  which  often  require  auxiliary  heaters  to  melt  the  sodium  to 
initiate  fluid  circulation. 

The  liquid  collects  in  the  bottom  of  the  evaporator  fins,  where  it 
contacts  the  wick  lining  the  inside  walls  of  the  fins.  The  wick  supplies 
liquid  sodium  over  the  entire  inside  surface  of  the  evaporator  through 
capillary  pumping. 

The  evaporator  configuration  was  based  on  the  more  difficult  heat 
transfer  from  the  flue  gas  to  the  evaporator  wall,  resulting  in  an  average 
heat  flux  of  100  kW/m2.  The  maximum  heat  flux  occurs  at  the  bottom  of  the 
evaporator  fins,  where  the  flue  gas  enters  the  heat  exchanger.  The 
magnitude  of  the  heat  flux  at  this  point  is  approximately  500  kW/mJ  .  In 
experiments  on  sodium  heat  pipes  at  Philips  and  STM,  this  magnitude  of  heat 
flux  was  defined  to  be  below  the  boiling  limitation  for  sodium  at  800°C, 
with  commonly  used  wire  mesh  wicks.  Once  the  configuration  was  established 
as  compatible  to  the  convective  and  radiative  properties  of  the  flue  gas  and 
the  boiling  properties  of  sodium,  the  design  task  was  directed  towards 
specification  of  a  wick  structure  capable  of  maintaining  an  adequate  supply 
of  liquid  sodium  over  the  entire  evaporator  surface,  i.e.  satisfying  the 
capillary  pumping  requirement,  and  providing  structural  support  to  the  flat 
fins. 

A  coarse  wire  mesh  was  placed  in  the  vapor  path  to  provide  structural 
support  for  the  fins  at  a  slight  penalty  in  temperature  uniformity  due  to 
vapor  flow  resistance,  while  simultaneously  providing  positive  contact 
between  the  capillary  wick  and  the  evaporator  wall.  A  schematic  cross- 
section  of  the  evaporator  is  shown  in  Figure  13. 


3.  Analysis  of  the  Heat  Exchanger/Evaporator 

An  analysis  was  conducted  on  the  evaporator  to  evaluate  capillary 
pressure  as  a  function  of  position  in  a  typical  fin  using  various  wick  and 
support  structures,  while  maintaining  the  external  geometry  defined  in  Table  2. 
The  analysis  was  developed  from  basic  heat  pipe  theory  and  empirical 
relations  for  flow  properties  of  wire  mesh  wicks  as  previously  summarized 
[10],  [11],  [12]. 

The  most  critical  aspect  of  the  evaporator  design  is  maintaining  suffi¬ 
ciently  high  capillary  pressure  in  the  wick  to  distribute  the  liquid 
throughout  the  fins.  In  addition,  the  coarse  supporting  mesh  must  withstand 
the  structural  load  due  to  the  pressure  difference  across  the  fin  walls 
without  significantly  restricting  the  vapor  flow  out  of  the  fin. 

Sonic  and  entrainment  limitations  are  a  much  lesser  concern  since  the 
large  number  of  fins,  hence  large  vapor  flow  area  results  in  vapor 
velocities  well  below  a  Mach  Number  of  0.2.  This  also  simplifies  the 
evaluation  of  the  vapor  pressure  drop  since  the  effects  of  the  vapor 
compressiblity  can  be  neglected. 

The  heat  flux  from  the  flue  gas  to  the  evaporator  varies  according  to 
the  temperature  distribution  as  shown  in  Figure  10.  However,  the  evaporator 
analysis  is  greatly  simplified  by  making  the  conservative  assumption  of 
constant  heat  flux.  This  assumption  results  in  liquid  flow  only  in  the  X 
(vertical)  direction,  and  vapor  flow  only  in  the  Y  (horizontal)  direction, 
essentially  reducing  the  analysis  to  a  one  dimensional  problem.  The  stream¬ 
lines  for  liquid  and  vapor  flow  in  both  cases  are  shown  in  Figure  14. 
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Figure  14a 

Fluid  Streamlines  in  a  Typical  Evaporator  Fin  with  Varying  Heat  Flux 
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Liquid  Streamlines 

Fin  Tip 


Fluid  Streamlines  in  a  Typical  Evaporator  Fin  with  Constant  Heat  Flux 


Application  of  Eq.  1  to  the  EHS  evaporator  results  in  the  following 


equation: 

AP  (x,y)  =  P  (ref)  +  AP  +  AP0  +  AP  (18) 

c  c  v  JC.  g 

=  capillary  pressure  at  position  (x,y) 

=  capillary  pressure  at  a  reference  position  (xrep  yre.) 
=  vapor  pressure  drop  between  (x,y)  and  (x  ,  y  ) 

=  liquid  pressure  drop  between  (x,y)  and  (x  f,  y  f ) 

=  hydrostatic  head  (due  to  gravitational  body  force) 

The  reference  position  was  defined  to  be  at  the  bottom  of  the  fins 
where  it  opens  into  the  evaporator  plenum  since  the  capillary  pressure  must 
be  at  its  minimum  value.  This  conclusion  can  be  deduced  readily,  since 
liquid  and  hydrostatic  pressures  here  are  defined  to  be  zero,  and  clearly 
the  vapor  pressure  must  be  at  its  smallest  level  in  the  fin  at  the  root. 

To  satisfy  Eq.  1  for  the  rest  of  the  heat  pipe,  the  minimum  capillary 
pressure  must  be  equal  to  the  sum  of  the  pressure  drops  in  the  rest  of  the 
heat  pipe. 

Vapor  and  liquid  flow  resistance  are  functions  of  the  mesh  geometry  and 
fluid  properties.  Development  of  the  liquid  and  vapor  pressure  drop 
equations  4  and  7,  respectively,  to  the  EHS  finned  evaporator  application 
can  be  found  in  Appendix  D.  The  results  of  the  analysis  are  summarized  in 
Figure  15. 

Several  simulations  were  conducted  using  various  meshes  fitting  the 
geometric  restrictions  as  noted  in  Table  2.  Several  wick  configurations 


where 


Pc(x,y) 

Pc(ref) 
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satisfied  the  capillary  pumping  requirement,  however,  the  chosen 
configuration  was  the  best  compromise  among  cost,  manufacturability  and 
performance.  The  results  of  the  performance  evaluation  are  summarized  in 
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Table  3. 


Table  3 


Heat  Exchanger/Evaporator  Construction  and  Performance  Data 


Number  of  fins 
Width/fin 
Wall  thickness 
Wick  structure 
Support  structure 


=  68 

=  2.67  mm 
=  0.2  mm 

=  4  layers  of  163  mesh,  .050  mm  (.002")  wire  diameter 
=  2  layers  of  8  mesh,  .48  mm  (.017")  wire  diameter 


Operating  Data 

Heat  pipe  evaporator  temperature  1085°K 

Heat  input  104  kW 

Maximum  wick  capillary  pressure  3200  Pa 

Performance 


Minimum  capillary  pressure  650  Pa 

Maximum  required  capillary  pressure  2525  Pa 

Temperature  drop  due  to: 

Conduction  across  evaporator  wall/wick  1.8  K 

Vapor  flow  restriction  in  fins  0.4  K 

Vapor  flow  restriction  in  connecting  tubes  2.0  K 

Condenser  (engine  heater)  4.8  K 

Total  temperature  drop  9.0  K 


Temperature  of  inside  of  wall  engine  heater  tubes  1076°K 


SECTION  IV 


EXTERNAL  HEATING  SYSTEM  CONSTRUCTION 

The  construction  of  the  EHS  was  designed  to  be  adaptable  to  simple 
folded  sheet  metal  manufacturing  techniques.  The  folded  constructions  are 
amenable  to  robotic  welding  techniques.  In  this  section,  a  brief 
description  of  the  envisioned  construction  procedure  is  presented. 

1.  Construction  of  the  EHS  Heat  Exchanger/Evaporator 

The  evaporator  construction  begins  by  bonding  the  layers  of  wire  mesh 
forming  the  wick  to  the  sheet  that  forms  the  outer  wall.  The  bonded  sheet 
is  folded  over  mandrels  of  1.0  mm  thickness  for  the  flue  gas  gap,  and  2.0  mm 
for  the  sodium  vapor  gap,  resulting  in  a  structure  like  that  shown  in  Figure  16. 
The  folded  sheet  is  then  soaked  in  solvent  to  remove  the  adhesive  and  any 
oil  or  grease. 

The  folded  sheet  is  then  mounted  on  a  special  mandrel  to  perform  the 
squaring  operation  on  the  root  of  the  fins  as  shown  in  Figure  17.  In 
addition,  the  edges  of  the  fins  are  crimped  together.  This  operation 
permits  ready  welding  of  the  fin  edges  without  filler  rod  and  simplifies  the 
welding  of  the  roots  of  the  fins  to  the  rest  of  the  enclosure. 

The  fins  are  filled  with  the  coarse  supporting  mesh,  then  mounted  in  a 
welding  fixture  with  heat  sinks.  At  this  stage,  the  fin  seams  are  welded. 

Once  the  fin  seams  are  welded,  the  assembly  is  placed  in  a  different 
welding  fixture  and  the  top  and  bottom  covers  are  welded  in  place.  These 
covers  fit  accurately  into  the  squared  fin  roots,  again  permitting  welding 
without  the  use  of  a  filler  rod.  See  Figure  18. 
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Figure  18 

Shape  and  Fit  of  Evaporator  Plenum  Walls  to  Fin  Roots 
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The  final  step  consists  of  filling  the  evaporator  plenum,  formed  by 
the  sides  of  the  outer  fins  and  the  top  and  bottom  covers,  with  coarse 
supporting  mesh,  then  welding  the  back  plate  in  place.  The  pack  plate 
contains  the  four  connecting  tubes,  with  mechanical  coupling  flanges,  that 
supply  the  STM4-120  ECU  with  sodium  vapor. 

2.  Construction  of  the  EHS  Recuperator 

The  techniques  used  to  form  the  recuperator  passages  are  identical  to 
the  initial  folding,  squaring  and  seam  welding  operations  used  to  form  the 
evaporator  enclosure. 

3.  EHS  Preliminary  Design  Drawings 

The  preliminary  design  drawings  for  the  complete  EHS  are  included  in 
Appendix  E. 

4.  Heat  Exchanger/Evaporator  Scale  Model  Tests 

In  order  to  confirm  the  functionality  of  the  novel  evaporator  design, 
scale  models  were  built  using  the  construction  techniques  outlined  above.  A 
photograph  of  the  model  is  shown  in  Figure  19. 


SECTION  V 


CONCLUSIONS 


The  conceptual  design  of  the  STM4-120  External  Heating  System  appears 
to  be  a  relatively  simple,  efficient,  and  inexpensive  system  for  supplying 
heat  to  the  Stirling  Energy  Conversion  Unit  to  provide  shaft  power  for  a 
mobile  electric  power  generator  set.  In  addition,  the  preliminary  design 
and  subsequent  scale  model  testing  of  the  EHS  heat  exchanger/evaporator 
substantiate  the  viability  of  the  conceptual  design. 

The  system  optimization  resulted  in  an  EHS  package  size  of  330  mm  x 
250  mm  x  70  mm,  operating  at  a  full  load  output  of  40  kW  with  an  equivalence 
ratio  of  1.25,  4%  atomizing  air,  fuel  consumption  rate  of  2.67  g/sec,  and  an 
EHS  exhaust  temperature  of  225°C. 

The  analysis  of  the  heat  exchanger /evaporator  and  scale  model  testing 
indicate  the  viability  of  the  evaporator  configuration  for  supplying  heat  at 
800°C  to  the  STM4-120  Stirling  engine.  The  analysis  also  provides  the 
framework  for  design  algorithms  for  future  designs. 

The  primary  concern  of  the  heat  exchanger/evaporator  preliminary  design 
was  to  specify  a  wick  structure  capable  of  providing  sufficient  liquid 
distribution  within  the  evaporator  fins  at  an  average  heat  flux  of  100  kW/m2 
and  a  maximum  of  500  kW/m2  . 

Results  of  the  preliminary  design  yielded  an  evaporator  with  68  fins, 
each  with  a  wick  structure  composed  of  four  layers  of  165  mesh,  .050  mm 
diameter  wire.  The  fin  support  structure  was  composed  of  two  layers  of  8 
mesh,  .43  mm  diameter  wire.  The  specified  wick,  capable  of  producing  3200 
Pa  of  capillary  pressure,  was  required  to  produce  only  2500  Pa  to  provide 
adequate  liquid  distribution. 


While  much  of  the  military  MEP  specifications,  MIL-STD-633  and 
MIL-G-52884  is  not  applicable  to  a  radically  new  and  different  prime  mover, 
such  as  the  STM4-120  Stirling  engine,  none  of  the  existing  design 
specifications  are  in  conflict  with  these  requirements. 

The  authors  feel  that  the  results  of  this  and  related  investigations 
warrant  further  development  of  the  STM4-120  for  Air  Force  MEP  applications. 
Potential  difficulties  not  accounted  for  in  the  preliminary  design  are 
problems  related  to  fouling  and  corrosion. 

Counter  flow  recuperative  preheaters  have  been  developed  by  various 
manufacturers  for  gas  turbines.  These  have  very  similar  requirements  to 
those  of  the  EHS  recuperator.  It  is  the  intention  of  STM  to  make  use  of  the 
know-how  of  the  gas  turbine  manufacturers,  particularly  in  the  areas  of 
fouling  and  oxidation  corrosion  of  the  recuperator  and  heat  exchanger/ 
evaporator,  and  also  in  addressing  weld  or  braze  joint  leakage. 

Furthermore,  it  is  possible  that  proven,  available  recuperators, 
metallic  or  even  ceramic,  can  be  directly  applicable  to  the  EHS  system. 

This  will  reduce  the  technical  risk  and  is  thus  the  preferred  approach 
provided  that  the  system  performance  does  not  have  to  be  excessively 
compromised. 

Solutions  to  these  problems  require  further  research  and  experimental 
testing  and  were  not  within  the  scope  of  this  effort. 

The  results  of  this  investigation  add  further  support  to  the 
conclusions  drawn  in  the  Argonne  study,  "Base  Technology  Stirling  Engine 
Military  Applications  Assessment,"  listed  below: 

•  Fuel  flexibility  (similar  to  that  possible  with  gas  turbine  sets); 


•  Low  specific  weight  (comparable  with  the  best  achievable  with 
diesel  sets); 

•  Low  noise  (better  than  that  achievable  with  silenced  diesel  sets); 

•  Low  IR  -  both  because  the  majority  of  heat  is  emitted  through  the 
radiator  at  already  low  temperature  and  because  the  engine  system 
is  more  efficient  than  current  sets  and  consequently  releases  less 
heat;  and 

•  High  efficiency  (comparable  to  the  best  achievable  by  diesels  in 
this  power  range. 

Stirling  Thermal  Motors,  Inc.  is  currently  engaged  in  the  fabrication 
of  five  prototypes  of  the  STM4-120  ECU  to  establish  the  predicted 
performance  and  reliability  of  this  Stirling  engine  configuration.  This 
program  is  expected  to  continue  into  1986,  during  which  prototypes  of  the 
EHS  design  are  anticipated  to  be  built  and  tested. 
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INTRODUCTION 


A  new  concept  in  Stirling  engine  technology  is  embodied  in  the 
engine  now  being  developed  at  Stirling  Thermal  Motors,  Inc.  in  Ann  Arbor, 
Michigan.  This  is  a  versatile  engine  suitable  for  many  different  applica¬ 
tions  and  heat  sources. 

The  engine,  rated  at  40  kW  at  2800  rpm,  is  a  four-cylinder,  double¬ 
acting  variable  displacement  Stirling  engine  with  pressurized  crankcase 
and  rotating  shaft  seal.  It  incorporates  indirect  heating  technology  with 
a  stacked  heat  exchanger  configuration  and  a  liquid  metal  heat  pipe 
connected  to  a  distinctly  separate  combustor  or  other  heat  source.  It 
specifically  emphasizes  high  efficiency  over  a  wide  range  of  operating 
conditions,  long  life,  low  manufacturing  cost  and  low  material  cost. 

This  article  describes  the  new  engine,  its  design  philosophy  and 
approach,  its  projected  performance,  and  some  of  its  more  attractive 
appl ications . 

BACKGROUND 

In  1072,  Ford  Motor  Company  obtained  a  worldwide  exclusive  license 
from  N.V.  Philips  of  the  Netherlands  for  the  Stirling  engine,  covering 
virtually  all  applications,  including  automotive. 

Under  this  license  agreement,  the  Research  Lab  at  N.V.  Philips 
was  to  design  and  build  four  175  hp  engines,  two  of  which  would  be 
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In  1976  the  two  Stirling  powered  Torinos  and  an  older  Philips 
Stirling  bus  equipped  with  a  4-cylinder  rhombic  drive  Stirling  engine  [2] 
were  successfully  demonstrated  for  three  days  in  Dearborn,  Michigan. 

A  few  years  later,  in  1978,  Ford  terminated  its  Stirling  engine 
activities  to  make  manpower  available  for  short-term  technological 
problems.  A  year  later.  Philips  also  stopped  work  on  the  Stirling  engine. 

Upon  these  events,  Stirling  Thermal  Motors,  Inc.  (STM)  was  founded 
in  the  United  States  to  continue  the  work  done  at  Philips,  so  that  the 
results  of  the  years  of  research  and  development  work  done  at  Philips 
since  its  last  license  agreement  in  1968  -  which  had  resulted  in  a 
technical  breakthrough  -  would  not  be  lost. 

STM's  main  purpose,  from  its  foundation,  was  to  develop  commercial 
Stirling  engines.  Philips  Laboratories  had  only  made  laboratory  models  for 
use  in  research.  The  only  engine  made  for  a  special  purpose  was  the  one 
for  Ford.  When  this  particular  engine  was  made.  Philips  was  confronted  with 
the  practical  reality  of  designing  and  building  a  Stirling  engine  for  the 
most  complex  application  imaginable  -  an  automotive  engine.  During  this 
time  it  was  discovered  that  some  components  of  the  engine  might  form 
obstacles  in  the  way  of  commercialization  because  of  their  complexity 
and  vulnerability. 

From  1974  on,  real  breakthroughs  were  made  in  avoiding  these 
complexities.  This  made  a  more  simple  four-cylinder,  double-acting  Stirling 
engine  possible. 

Unfortunately,  by  this  time  it  was  too  late  to  incorporate  these 


improvements  into  the  Ford  engine.  The  intent  was  to  use  these  new 
developments  in  a  second-generation  Ford  engine.  This,  though,  was  not 


done  before  Ford  dropped  its  Stirling  engine  program.  The  engines  being 
built  by  earlier  Philips  licensees  were  based  on  designs  older  than  the 
Ford  engine.  Their  configurations  had  been  frozen  for  several  years.  It 
was  therefore  impossible  to  utilize  the  new  improvements. 

From  the  outset  STM  was  convinced  that  the  time  was  ripe  for 
commercialization  of  the  Stirling  engine  because  all  the  ingredients 
for  a  simple,  inexpensive  and  reliable  engine  with  a  long  service  life 
were  present. 


GENERAL  APPROACH,  BASIC  APPLICATIONS 


STM's  general  approach  is  based  on  the  conclusion  that  competition 
with  existing  internal  combustion  engines  should  be  avoided,  at  least  in 
the  beginning.  Rather,  markets  should  be  found  where  the  IC  engine  cannot 
be  used  and  where  the  use  of  the  Stirling  engine  would  be  very  economical, 
making  use  of  the  unique  properties  of  the  Stirling  engine.  Of  the  many 
possible  applications,  particular  attention  was  given  to  the  following 
three: 


•  Prime  mover  for  heat-driven  heat  pump 

•  Solar  energy  conversion 

•  Engine  for  generator  sets. 


If  the  manufacturing  cost  of  the  engine  could  be  sufficiently  low, 
particularly  in  mass  production,  the  market  in  these  fields  alone  could  be 


TECHNICAL  APPROACH 


The  whole  drive  of  STM  is  to  commercialize  the  Stirling  engine.  This 
means  that  the  engine  must  be  simple,  reliable,  inexpensive,  and  that  it 
should  have  a  long  service  life.  None  of  these  requirements  should  have  an 
adverse  effect  on  the  performance  of  the  engine. 

More  than  three  years  of  designing,  discussions  with  suppliers  and 
vendors,  component  testing  and  price  calculations,  led  to  the  new  version 
of  the  Stirling  engine.  Studies  done  for  NASA  have  shown  that  this  engine 
configuration  is  suitable  for  a  whole  range  of  power  sizes  up  to  500  hp. 

Special  emphasis  was  placed  on  the  flexibility  of  the  engine  to 
adapt  readily  to  a  wide  range  of  specific  applications,  duty  cycles  and 
heat  sources. 

Consequently,  indirect  heating  technology  is  an  integral  part  of 
the  development  effort,  making  it  possible  to  divide  the  engine  into  an 
energy  conversion  unit  (ECU)  and  a  distinctly  separate  external  heating 
system  (EHS).  Different  heat  sources  coupled  to  the  same  ECU  will  adapt 
the  engine  to  different  applications  and  enhance  commercial  introduction 
since  most  of  the  development  complexity  and  cost  is  in  the  energy  con¬ 
version  unit. 

A  liquid  metal  heat  pipe  [3]  is  used  to  transport  the  heat  from  the 
heat  source  to  the  expansion  heat  exchanger  of  the  thermal  converter. 

So  far,  most  of  the  development  effort  has  concentrated  on  the  ECU, 
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which  is  designated  STM4-120  (4  cylinders,  120  cm  swept  volume  per  piston) 
This  engine  is  shown  in  Figures  3a  and  3b. 
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ADVANTAGES  OF  REMOTE  HEATING 


One  of  the  obstacles  in  the  way  of  mass  production  of  the  Ford/Philips 
engine  was  the  heater  head.  This  was  built  as  an  integrated  unit  for  the 
four  cylinders  (Figure  4).  The  huge  mass  of  heat  resistant  material  was  very 
expensive  and  made  the  brazing  cycle  much  too  long.  The  reason  for  this 
large  amount  of  material  is  that  the  tubular-expansion  heat  exchanger  common 
to  direct  flame  Stirling  engines  must  accomodate  the  relatively  difficult 
heat  transfer  from  the  flue  gas  to  the  walls  of  the  heat  exchanger  tubes. 


Figure  4 

The  integrated  direct-f la«e-heated  heater  head  (fro«  the  Ford/Philips  engine). 


It  is,  therefore,  characterized  by  a  complex  cage  geometry  as  well  as  volume 
and  flow-path  length  which  are  much  larger  than  those  required  for  the 
relatively  easy  heat  transfer  from  the  tube  walls  to  the  working  fluid 


By  contrast,  an  expansion  heat  exchanger  heated  by  the  condensing 
metal  vapor  with  a  large  film  coefficient  in  a  heat  pipe  can  be  ideally 
sized  to  suit  the  requirements  of  the  working  fluid  and  can  be  shaped  in 
the  most  convenient  manner  for  ease  of  fabrication  (Figure  5). 


Figure  5 

Expansion  heat  exchangers  of  the  STM4-120.  There  is  one  per  cylinder.  These  will  later  be  electron-beam 
welded  in  the  heat  exchanger  stack.  The  tubes  are  curved,  enclosed  in  a  flexible  cannister,  and 
brazed  to  two  end  plates. 


Of  course,  this  itself  does  not  solve  the  difficult  external  heat 
transfer  problem,  but  rather  shifts  it  to  the  evaporator  section  of  the  heat 
pipe  where  the  size  necessary  for  adequate  heat  transfer  is  easily 
realized  since  the  heat  pipe  does  not  have  to  support  the  high  cycle 


pressure. 


Indirect  heating  thus  offers  a  number  of  advantages  in  addition  to 


the  flexibility  with  which  it  endows  the  engine: 

•  It  brings  about  major  simplifications  to  the  heat-exchanger 
design.  The  so-called  heat-exchanger-stack  configuration, 
designed  to  take  advantage  of  the  high  film  coefficient  of  the 
condensing  metal  vapor,  is  considerably  less  expensive  and 
more  suitable  for  mass  production. 

•  It  brings  about  considerable  improvement  of  the  engine 
performance  by  permitting  the  heater  design  to  be  ideally 
suited  to  the  thermodynamic  requirements. 

•  The  uniform  temperature  throughout  the  confines  of  the  heat 
pipe  enclosure  eliminates  hot  spots  on  the  heater  and 

thus  enhances  both  the  efficiency  and  the  reliability  of 
the  engine. 

NEW  POWER  CONTROL  SYSTEM 

Up  to  this  time,  the  preferred  method  for  changing  power  was 
changing  the  pressure  inside  the  engine,  because  the  torque  of  the  engine 
is  approximately  proportional  to  the  mean  pressure  of  the  working  gas  [4], 
The  development  of  this  type  of  power  control  at  Philips  was  done 
with  a  single  cylinder  displacer  engine,  where  this  type  of  power  control 
was  acceptable.  However,  for  a  four-cylinder,  double-acting  engine  it 
became  quite  cumbersome,  particularly  when  very  rapid  changes  were  required, 
as  in  automotive  applications.  This  type  of  system  included  many  check 


valves,  activator  valves,  and  a  storage  bottle,  along  with  a  high  pressure 
hydrogen  compressor. 

Figure  6  shows  a  diagram  of  the  power  control  system  of  the 
Ford/Philips  Stirling  engine.  Power  increases  when  the  working  gas  (in 
this  case  hydrogen)  is  dumped  from  the  high  pressure  storage  bottle  into 
the  engine.  The  reverse  takes  place  when  the  gas  is  pumped  out  of  the  engine 
into  the  storage  tank  with  the  high  pressure  compressor.  But  because  this 
is  a  slow  process,  during  this  time,  a  short-circuit  power  control  -  which  is 
a  loss  control  -  instantaneously  cuts  the  power. 
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Figure  6 

Power  control  system  for  the  Ford/Philips  Stirling  engine,  the  torque  of  the  engine  is  controlled  by 
the  pressure  of  the  working  gas.  For  more  power,  the  working  gas  from  the  storage  bottle  is  supplied  t 
the  engine.  For  less  power,  small  hydrogen  compressors  (connected  to  the  bottoms  of  the  pistons)  are 
pumping  the  gas  out  of  the  engine  and  back  into  the  storage  bottle. 


In  1974,  during  the  work  on  the  automotive  engine,  a  relatively  simple 


heavy  duty  construction  was  found  to  vary  the  power  [5],  [6],  [7],  In  this 
case  the  mean  pressure  of  the  engine  stays  the  same,  but  the  stroke  of 
the  pistons  changes.  This  method  of  power  control  has  the  further  advan¬ 
tage  of  high  part  load  efficiency.  Such  a  construction  could  be  used  only 
with  a  swashplate  drive  since  the  stroke  of  the  pistons  is  controlled  by 
the  angle  of  the  swashplate.  It  was  tested  thoroughly  in  a  test  rig  and 
applied  in  the  Advenco  engine,  but  the  Advenco  engine  was  never  adequately 
tested.  Philips  eventually  sold  it  to  NASA,  where  further  testing  was  done. 

Figure  7  illustrates  schematically  the  variable  swashplate  mechanism. 
The  swashplate  is  mounted  on  a  part  of  the  shaft  which  is  tilted  an  angle  a 
from  the  main  shaft  axis.  The  swashplate  is  mounted  in  such  fashion  that  its 
centerline  makes  an  angle  a  with  the  tilted  shaft  axis  and  it  can  be  rotated 
relative  to  and  about  the  tilted  shaft  axis  in  order  to  change  its  angle  and 
with  this,  the  stroke  of  the  pistons. 


figure  7  -  Schematic  of  the  variable  swashplate  mechanism,  showing  the  principle  of  changing  the  angle  of 
the  swashplate,  which  It  makes  with  the  shaft  from  0  to  2  a.  The  small  plate  A  Is  fixed  on  the  shaft  S  with 
an  angle  a  .  The  engine  swashplate  (drawn  here  as  a  solid  line)  Is  perpendicular  on  the  shaft  S.  This  situation, 
a  -  n  »  0,  means  that  the  strokes  of  the  pistons  are  zero.  When  the  engine  swashplate  B  Is  turned  180®  relative 
to  the  small  plate  A,  the  angle  then  becomes  represented  by  the  dotted  lines.  In  this  case,  the  strokes  of  the 
pistons  are  maximal.  By  turning  8  relative  to  A  between  0®  and  180®  any  angle  of  the  swashplate  between  0  and  2a 

can  be  obtained,  so  the  strokes  of  the  pistons  can  be  changed  from  zero  to  maximum.  C  Is  a  hydraulic  vane 

mechanism,  the  housing  of  which  can  turn  In  one  or  the  other  direction  depending  on  which  chambers  are  pressur¬ 
ized  and  which  are  not.  The  oil  is  fed  via  two  channels  In  the  shaft  S.  The  turning  of  the  housing  Is  trans¬ 

mitted  to  the  engine  via  bevel  gears  E  and  F. 


Figure  9 

Practical  model  of  a  variable  swashplate,  shown  in  two  positions. 

The  torque  applied  by  the  actuator  to  the  swashplate  in  order  to 
maintain  a  certain  angular  position  depends  on  whether  it  was  rotated 
to  such  position  in  the  positive  or  in  the  negative  direction.  Rotation 
in  the  negative  direction  requires  less  torque  since  the  engine  torque 
itself  acts  in  this  case  to  increase  the  swashplate  angle. 

Figure  10  shows  the  actuator  torque  as  a  function  of  the  swashplate 
angular  position  relative  to  the  shaft.  The  curves  labeled  M+(4>)  and 
M  (ij;)  refer  to  that  torque  for  the  positive  and  negative  direction  of 
rotation  respectively.  The  third  curve,  y(^)  shows  the  corresponding 
swashplate  angle.  The  curve  reverses  its  sign  within  its  range  of 

definition.  The  point  of  sign  reversal  is  an  unstable  control  point  to 
be  avoided  by  narrowing  the  range  of  definition  to  exclude  it.  In  the 


Figure  10 

Rotary  actuator  torque  requirements.  H+  means  relative  rotation  for  larger  stroke  in  the  same  direction 
as  the  rotation  of  the  engine.  means  relative  rotation  for  larger  stroke  in  the  opposite  direction 
from  the  rotation  of  the  engine. 


Energy  Conversion  Unit,  the  angle  a  is  12.5°,  yielding  maximum  theoretical 
swashplate  angle  of  2a  =  25°,  corresponding  to  180°  rotation.  The  maximum 
swashplate  angle  of  interest  is  only  22°,  corresponding  to  a  narrower  range 
of  definition  (124°)  in  which  the  torque  M  (ij/)  does  not  change  sign. 

Loss  of  hydraulic  power  will  result  in  the  gas  forces  bringing  the 


swashplate  to  a  position  perpendicular  to  the  main  shaft  axis,  reducing  the 
piston  stroke  to  zero  -  an  automatic  safety  feature. 


In  a  four-cylinder,  double-acting  Stirling  engine  there  are  two 


types  of  dynamic  seals: 

•  Dynamic  seals  as  piston  rings  to  divide  the  four  cycles  from 
each  other,  and 

•  Dynamic  reciprocating  seals  on  the  piston  rods,  in  order  to 
contain  the  high  pressure  working  gas  in  the  engine.  These 
seals  should  also  prevent  oil  penetration  into  the  cylinders 
from  the  lubricated  drive. 

For  the  dry-running  piscon  rings,  a  good  solution  is  found  in 
using  a  reinforced  PTFE  material. 

However,  the  different  types  of  reciprocating  seals  for  the  piston 
rods  are  still  not  reliable.  Philips  developed  the  rolling  diaphragm  seal, 
but  this  was  shown,  in  the  Ford/Philips  engine,  to  be  vulnerable  in 
non-laboratory  environments. 

STM  was  able  to  avoid  the  gas  containment  function  of  these 
reciprocating  seals  entirely. 

The  new  power  control,  with  its  variable  swashplate,  made  it 
possible  to  enclose  the  relatively  small  drive  with  a  pressure  hull  and  to 
use  a  commercially  available  rotating  shaft  seal.  Preventing  oil  penetration 
into  the  cylinders  is,  in  this  case,  much  easier  and  has  already  been 
thoroughly  tested  in  other  engines. 
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SPECIAL  FEATURES 


Amenability  to  dynamic  balancing  and  the  ease  of  starting  the 
engine  are  two  additional  features  of  the  variable  swashplate  drive 
and  power  control  elaborated  upon  in  this  section. 

DYNAMIC  BALANCING  is  achieved  by  adjusting  the  swashplate  moments 
of  inertia  to  the  reciprocating  mass.  This  is  done  in  a  manner  enforcing 
perfect  dynamic  balance  at  a  certain  swashplate  angle  within  its  range 
of  variation.  At  different  angles  unbalanced  moments  will  appear,  but  since 
perfect  balance  automatically  occurs  at  zero  angle,  these  will  be  very 
smal 1  . 

STARTING  of  the  engine  can  be  accomplished  by  heating  up  the 
expansion  heat  exchanger  and  the  regenerator  and  then  suddenly  using  the 
control  pressure  to  increase  the  swashplate  angle.  This  will  cause 
the  pistons  to  move  in  their  normal  way,  causing  the  engine  to 
immediately  develop  sufficient  power  to  perpetuate  the  motion.  An 
accumulator  fed  by  the  hydraulic  pressure  pump  will  be  used  for 
that  purpose. 

Obviously,  such  a  simple  procedure  may  be  used  only  for  such 
applications  as  solar  conversion  since  no  accessories  are  required  for  the 
combustion.  In  other  cases  only  a  very  small  starter  motor  is  required 
to  power  the  accessories  needed  for  combustion,  such  as  the  air  blower. 

When  the  engine  has  reached  its  correct  temperature  the  accumulator 
pressure  may  be  used  to  quickly  increase  the  angle,  having  the  engine 
self-start  as  mentioned  above. 


DESCRIPTION  OF  THE  ENGINE 


A  layout  drawing  of  the  Energy  Conversion  Unit  is  shown  in  Figure  11. 
This  unit  is  distinguished  by  two  major  features: 

•  Variable  swashplate  drive  and  power  control  contained  in  a 
pressurized  "crankcase"  with  a  commercially  available  rotating 
shaft  seal  containing  the  working  fluid  and  making  it  possible 
to  avoid  the  reciprocating  rod  sealing  problem;  and 

•  Indirect  heating,  featuring  heat-exchanger-stack  configuration 
and  liquid  metal  heat  pipe  heat  transport  system. 

Following  is  a  brief  description  of  key  components. 

The  rotating  shaft  seal  assembly  (1)  is  bolted  to  the  rear  crankcase. 
The  drive  shaft  (2)  is  splined  to  the  main  shaft  to  allow  for  slight 
deflection  or  misalignment  without  affecting  the  seal  performance.  The 
drive  shaft  is  supported  radially  by  one  needle  and  one  ball  bearing,  but 
since  it  passes  through  the  pressurized  crankcase,  it  must  support 
the  corresponding  thrust  load.  The  thrust  bearing  (3)  is  a  heavy  duty  tapered 
roller  bearing. 

The  rotating  shaft  seal  (4)  is  a  commercially  available  mechanical 
face  seal  similar  to  those  found  in  compressors  and  gearboxes.  The  seal  is 
secured  in  the  housing  by  an  interference  fit.  A  rotating  mating  ring  (5) 
is  mounted  on  the  shaft.  Both  components  have  an  extremely  flat,  smooth 
sealing  face.  When  the  seal  is  installed,  these  faces  are  pressed  together 
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Figure  11 

The  STM4-120 

1) 

Rotating  shaft  seal  assembly 

10)  Rotary  actuator 

19) 

Heater 

2) 

Orive  shaft 

11)  Main  shaft 

20) 

Hot  connecting  duct 

3) 

Rotating  shaft  seal 

12)  Bevel  gear 

21) 

Piston  assembly 

A) 

Mating  ring 

13)  Swashplate 

22) 

Cylinder/regenerator  housing 

5) 

Radial  lip  seal 

14)  Front  main  bearing 

23) 

Hydraulic  service  assembly 

6) 

Thrust  bearing 

15)  Oil  pump  module 

24) 

Power  control  valve 

7) 

Pressure  hull 

16)  Front  crankcase 

25) 

Oil  scraper/capseal  assembly 

8) 

Rear  main  bearing 

17)  Cooler 

26) 

Crosshead 

9) 

Rear  crankcase 

18)  Regenerator 

27) 

Accumulator 

and  held  in  contact  by  compressed  springs  in  the  seal,  then  surrounded  by 
oil.  Because  of  their  closely  controlled  surface  finish,  virtually  no 
oil  can  pass  between  the  sealing  faces. 

To  minimize  wear,  the  seal  is  hydraulically  balanced  to  reduce  contact 
forces  due  to  crankcase  pressure  and  subjected  to  a  steady  flow  of  oil 
to  maintain  cool  operation. 

The  radial  lip  seal  (5)  prevents  the  oil  surrounding  the  rotating 
shaft  seal  from  draining  into  the  sump. 

The  pressure  hull  (7)  fits  over  the  rotating  shaft  seal  housing 
and  is  bolted  to  the  front  of  the  crankcase,  covering  the  entire  drive. 

The  pressure  hull  serves  no  structural  function  other  than  containing  the 
crankcase  pressure. 

The  front  (16)  and  rear  (9)  crankcase  castings  form  the  major 
building  blocks  of  the  engine,  housing  the  entire  drive  and  power  control. 

In  addition,  the  front  crankcase  houses  the  cold  thermodynamic  system 
including  the  coolers,  water  passages  and  cold  connecting  duct.  Both 
are  cast  from  Reynolds  390  Aluminum  for  reduced  weight  and  suitable  running 
surfaces  for  the  reciprocating  components. 

The  swashplate  (13),  combined  with  the  main  shaft  (11)  provide 
the  mechanism  for  converting  the  reciprocating  motion  of  the  pistons  into 
rotary  motion.  The  swashplate  is  mounted  on  the  shaft  in  a  manner  such 
that  relative  rotation  to  the  shaft  causes  the  swashplate  angle  to  change. 

The  front  (14)  and  rear  (8)  main  bearings  provide  support  for  the 


The  rotary  actuator  (10)  is  a  hydraulic  vane  motor  comprising  two 
diametrically  opposed  vanes  attached  to  the  shaft  and  two  attached  to  the 
housing,  forming  two  chambers.  Rotation  of  the  actuator  housing  relative 
to  the  shaft  is  accomplished  by  pressurizing  one  chamber  and  relieving  the 
other.  The  relative  rotation  between  the  shaft  and  actuator  housing  is 
transmitted  to  the  swashplate  via  the  bevel  gear  (12)  to  change  the 
swashplate  angle  and  therefore  the  piston  stroke. 

The  oil  pump  module  (15)  houses  the  lubrication  and  hydraulic  pumps 
which  are  directly  driven  by  the  main  shaft.  Both  pumps  are  high-efficiency 
geroters.  Internal  passages  duct  the  oil  to  galleries  to  provide  lubrication 
High  pressure  oil  is  ducted  to  the  externally  mounted  hydraulic  control 
system  and  through  annular  passages  in  the  shaft  to  the  rotary  actuator. 

The  hydraulic  service  assembly  (23)  is  mounted  to  the  underside  of 
the  front  crankcase  for  ready  access.  It  consists  of  the  oil  sump,  hydraulic 
accumulator  and  the  power  control  valve. 

The  accumulator  (27)  maintains  a  large  volume  of  oil  at  hydraulic 
system  pressure  to  prevent  lags  in  response  to  power  demands.  The  accumulator 
permits  use  of  a  smaller  pump,  thus  reducing  parasitic  losses. 

The  power  control  valve  (24)  is  a  solenoid  activated  proportional 
spool  valve  that  controls  the  pressure  difference  across  the  rotary  actuator 
and  in  turn,  the  angle  of  the  swashplate  and  piston  stroke. 

The  four  coolers  (17)  are  conventional  tubular  heat  exchangers  mounted 
in  the  front  crankcase.  There  is  one  cooler  for  each  cycle. 

The  high  temperature  components  consist  of  identical  assemblies  for 
each  cycle  forming  a  heat  exchanger  stack  configuration.  These  components 
are  investment  cast  from  the  lower  cost  iron-based  superalloy  CRM-6D. 


The  engine  heater  (19)  is  a  simple  tube  bundle  heat  exchanger  consisting 
of  relatively  short  and  small  diameter  tubes  enclosed  in  a  flexible  cannister 
and  brazed  to  end  plates.  The  heater  is  then  electron  beam  welded  to  the 
hot  connecting  duct  (20)  and  cyl inder/regenerator  housing  (22). 

After  installation  of  the  regenera  tor  (18),  a  sintered  stack  of  stain¬ 
less  steel  wire  mesh,  each  assembly  is  mounted  on  the  front  crankcase. 

The  piston  (21)  is  a  light  weight  welded  assembly  consisting  of  a  base, 
which  houses  the  dry  running  piston  rings  and  bearing,  welded  to  the  piston 
rod  and  dome. 

The  oil  scraper/capseal  assembly  (25)  separates  each  cycle  from  the 
crankcase.  The  capseal  separates  the  fluctuating  cycle  pressure  from  the 
mean  crankcase  pressure.  The  oil  scraper  is  a  babbit  sleeve  radially  pressed 
against  the  rod.  There  is  no  pressure  difference  across  the  scraper.  Its 
only  function  is  to  prevent  oil  from  penetrating  into  the  cycle.  Both  the 
oil  scraper  and  the  capseal  are  mounted  in  a  compliant  housing  to 
accomodate  the  slightly  eccentric  motion  of  the  piston  rod  due  to  the 
running  clearance  of  the  cross  head. 

The  cross  heads  (26)  are  made  from  a  high  strength  aluminum  alloy. 

The  extended  length  reduces  side  forces  and  bridge  deflection.  The  running 
surfaces  are  treated  with  a  material  that  is  compatible  with  the  Reynolds 
390  bore. 

The  life  of  the  engine  is  determined  by  the  creep  failure  life  of  the 
hot  parts,  and  is  well  in  excess  of  50,000  hours  at  full  load. 

Table  1  summarizes  some  of  the  important  features  and  parameters 
of  the  engine. 
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Arrangement: 

Bore: 

Maximum  Stroke: 

Overall  Length: 
Cross-sectional  Dimensions: 
Total  Estimated  Weight: 
Working  Fluid: 

Mean  Cycle  Pressure: 

Heater  Temperature: 

Power  Control : 

Teat  Transport: 

Gas  Containment: 

Oil  Containment: 

Material s : 


Four  double-acting  cylinders  symmetrical ly 
arranged  about  a  common  axis.  One  heat 
exchanger  assembly  per  cylinder. 

56  mm 

48  mm 

635  mm 

Largest  cross  section  is  300  mm  in  diameter 

75  kg 

Hel ium 

11  MPa 

800°C 

Piston  stroke  variation  by  means  of  a 
variable  swashplate  with  a  maximum  angle 
of  22° 

Sodium  heat  pipe 

Crankcase  pressurized  to  mean  cycle  pressure 
and  sealed  with  a  rotating  shaft  seal 

Reciprocating  rod  oil  scrapers 

Iron-base  CRM-6D,  CG-27  heater  tubes 


Table  1 


Important  Features  and  Parameters  of  the  STM4-120 


PERFORMANCE 


A  set  of  engine  parameters  governing  the  performance  of  the 
thermodynamic  section  was  selected  for  the  STM4-120  through  an  extensive 
and  painstaking  effort  to  tailor  the  engine  performance  to  the  design 
approach  described  above. 

This  approach  required  a  high  level  of  efficiency  to  prevail  over 
a  wide  range  of  operating  conditions,  making  the  engine  suitable  for 
any  duty  cycle.  Hundreds  of  different  combinations  of  engine  parameters 
were  simulated  before  the  optimal  set  was  selected  and  established  as 
the  design  base  of  the  engine. 

The  result  is  shown,  in  Figures  12a  and  12b,  as  two  performance  maps 
at  mean  cycle  pressure  levels  of  11  MPa  and  6.3  MPa  respectively.  Figure  12a 
shows  that  between  the  power  levels  of  8  kW  and  40  kW  the  shaft  efficiency 
is  between  45%  and  47%  (excluding  auxiliaries  and  collector  efficiency). 

Over  a  wide  power  range  the  efficiency  hardly  changes  with  engine  speed. 

For  applications  requiring  less  power,  a  smaller  charge  of 
helium  may  be  used  with  very  little  effect  on  the  efficiency.  This  is 
shown  in  Figure  12b  for  an  engine  charged  to  6.3  MPa  to  provide  power 
output  of  no  more  than  25  kW.  This  will  greatly  enhance  the  life  of  the 
engine. 

Major  contributions  to  the  high  efficiency  of  the  engine  are  derived 
from  the  heat-exchanger-stack  configuration,  the  variable  stroke  power 
control  and  the  design  of  the  drive  mechanism. 

The  heat-exchanger-stack  configuration  adds  more  freedom  to  the  set 


of  parameters  governing  the  thermodynamic  performance.  These  can  be 
exploited  to  facilitate  tailoring  of  the  performance  characteristics  of 
the  engine. 


The  variable  stroke  power  control  inherently  inhibits  degradation 
of  the  efficiency  at  part  load  since  power  reduction  is  accomplished  partly 
through  the  addition  of  void  volume  or,  equivalently,  reduction  of  the 
pressure  wave  amplitude,  which  is  beneficial  to  the  indicated  efficiency. 

Figure  13  schematically  illustrates  the  ECU  with  solar  receiver, 
solid  fuel  combustion,  and  with  a  liquid  or  gaseous  fuel  combustor. 


Solar  Receiver 


Figure  13 

The  STH4-120  with  different  heat  sources 
80 
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Symbols 


Heat  transfer  area  [m2]  (for  gap:  dA  =  2wdx) 

Gap  width  [m] 

Heat  transfer  coefficient  [  ] 

m2  °C 

Thermal  conductivity  [  ^  ] 

Mass  flow  [kg/sec]  per  gap 
Nusselt  number  (for  gap:  Nu  =  hg/k) 

Heat  flow  (Watts) 

Temperature  (°C) 

Preheater  plate  wall  thickness  (m) 

Depth  of  heat  exchanger  (or  preheater)  gap  (m) 

Coordinate  in  the  flue  gas  flow  direction  (m) 

Specific  enthalpy  (3/kg) 

Subscripts 

f  denotes  properties  of  the  flue  gas 

a  denotes  properties  of  the  incoming  air 

w  denotes  properties  of  the  heat  exchanger  wall  material 


The  analysis  is  performed  on  a  counterflow,  flat  plate  recuperator  with 
multiple  parallel  gaps  of  width  (g),  depth  (w)  and  length  (  It  ). 

Laminar  flow  is  assumed  throughout  the  gaps  and  hence  the  Nusselt 
number  (based  on  the  gap  width)  is  constant: 


Nu  =  ^  g  =  3.8  (1) 

The  analysis,  performed  on  a  single  gap,  takes  into  account  the  varia¬ 
tion  of  fluid  properties  with  temperature,  but  not  with  pressure. 

The  functional  relationship  of  the  thermal  conductivity  (k)  and  the 
enthalpy  (T)  to  the  temperature  T  for  both  flue  gas  and  incoming  air  is 
known  for  the  particular  fuel  used  and  the  equivalence  ratio.  The  latter 
relation  takes  into  account  molecular  dissociation. 

The  energy  balance  for  the  flue  gas  requires  that  heat  loss  from  the 
flue  gas  by  convection  and  conduction  to  the  incoming  air  stream  be  equal  to 
the  reduction  in  enthalpy  of  the  flue  gas  stream: 
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k  (T  )  1  +k 
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d  cpf(Tf) 
dx 


+  2w(Tf  -  Ta)  =  0  (3) 
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Similarly,  energy  balance  for  the  air  stream  requires  that: 


dQ  =  m  d  <p  (T  ) 

d  dad 


but 


hence: 


dQa  =  -dQf  =  mfd<pf(Tf) 


W 


(5) 


•  '“Ca'V 

m  - — a - 

a  dx 


=  -m 


.  dcpf(T{) 


f  dx 


(6) 


Equations  (3)  and  (6)  form  a  set  of  two  coupled,  first  order,  non¬ 
linear  ODE's  for  the  two  unknown  functions  T,(x)  and  T  (x).  The  boundar 

f  a 

conditions  are  the  known  inlet  temperatures  T^(O)  and  Ta(£)  of  the  flue  gas 


and  air  streams  respectively. 

Governing  Equations  for  the  Heat  Exchanger/Evaporator  Analysis 

The  analysis  is  performed  on  a  heat  exchanger  with  flat  parallel  gaps 

(of  width  g,  depth  w,  and  length  1  )  through  which  flue  gas  flows  laminarly 

at  a  constant  Nusselt  number  Nu  --  3.8.  The  gas  yields  its  heat  to  the  gap 

walls  which  are  at  constant  and  uniform  temperature  t  . 

w 


V  V  -v. 
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As  before,  energy  balance  on  the  flue  gas  requires: 


Substitution  of 


dQf  =  -hf(Tf  -  Tw)  dA  -  mfdcpf(T{) 


and  of 


in  (7)  yields: 


dA  =  2wdx 


.  dcp(T)  k(T) 

m.  +  2w  Nu—  -C  (T,  -  T  )  =  0 

f  dx  g  f  w 


This  first  order,  non-linear  ODE  in  the  unknown  function  T^(x)  forms  the 
mathematical  model  of  the  evaporator/heat  exchanger.  The  boundary  condition 
is  the  known  inlet  temperature  Tf(o)  which  is  the  flame  temperature. 


p 


Fuel  Data 

HTCR 

SLHV 


EHS  Analysis  Code  Listing 

atomic  ratio  of  hydrogen  to  carbon  in  fuel 
specific  lower  heat  value  of  fuel,  cal./gr. 


Operating  Data 

ER  equivalence  ratio 

AAR  atomizing  air  to  intake  air  ratio 

TPA  initial  guess  for  preheated  air  temperature, 

MFUEL  fuel  consumption  rate,  gr./sec. 

TWEX  outside  wall  temperature  of  evaporator  °K 

TAIN  air  intake  temperature  °K 


Geometric  Data;  Heat  Exchanger  (see  figure) 


EXL 

length,  m 

WEX 

width,  m 

NGEX 

number  of  flue  gas  passages  ("gaps 

") 

GEX 

gap  width,  m 

NUEX 

Nussel  number 

THM 

wall  thickness,  m 

CM 

thermal  conductivity  of  the  metal. 

w/m  °K 

Geometric 

Data;  Preheater  (see  figure) 

PL 

length,  m 

WP 

width,  m 

NGP 

GPA 

GPF 

NUP 


number  of  flow  passages  for  each  stream 
gap  width,  air  side,  m 
gap  width,  flue  side,  m 
Nusselt  number 


Cold  Start  Data 


ROM 

CAPM 

ETHG 


density  of  the  metal,  gr/cm^ 

specific  heat  of  the  metal,  cal./gr.  °K 

effective  gauze  thickness:  the  thickness  of  a  solid  sheet 
of  metal  of  the  same  weight  per  unit  area  as  the  gauze 
stack,  m 
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number  of  discrete  steps  in  heat  exchanger  calculation 

number  of  discrete  steps  in  preheater  calculation 

maximum  convergence  error  in  preheater  iterative 
procedure,  relative 

maximum  convergence  error  in  EHS  iterative  procedure 
relative 


SAMPLE  RUNS  -  AUGUST  8,  1985 


1.8040  10179. 

1.2500  . 40000E-01  949.90  2.0000  1075.0  301.00 

.12000  . 60000E— 0 1  61  . 13000E-02  3.8000  . 30000E-03  21.600 

.18000  . 60000E-01  78  . 13000E-02  . 13000E-02  3.8000 

8.0000  .12000  .  00000, 

100  100  . 20000E-02  . 20000E-02 

PREHEATED  AIR  TEMP.  949.9  K 

FLAME  TEMPERATURE  2374.9  K 

HEAT  EXCHANGER  OUTPUT  TEMP.  1096.3  K 

HEAT  OUTPUT  74.661  kw 

EXHAUST  OUTLET  TEMP.  569.9  K 

E.H.S.  EFFICIENCY  87.675 

COLD  START  PENALTY  45.389  gr 

PRESSURE  DROP  2.193  K Pa 

FLOW  FRICTION  POWER  296.791  w 

623.38  841.16  1388.4 


B>A: PRINT  EHS.FOR 


Nat  ready  error  reading  drive  A 
Abort,  Retry,  Ignore?  R 


BsEHS 


.FOR  is  currently  being  printed 


3s  0.12,0.06,  61  1 . 3E-3, 3. 8, O. 3E-3, 21 . 6, 

4:  0.18,0.06,  78  1 . 3e-3, 1 . 3e-3, 3. 8, 

5s*  8.0,0. 12,0. 537E-3, 

6:  100  100  0.002,0.002 


B>EHS 
1.  .  8 04 A 
1 . 2 50 A 
. 12000 
. 18000 
8. 0000 


10179. 

. 40000E-01  949.90  2.0000  1075.0  301.00 

. 60000E-01  61  . 13000E-02  3.8000  . 30000E-03  21.600 

. 60000E-01  78  . 13000E-02  . 13000E-02  3.8000 

. 1 2000  . 5370QE-03 


.18000  . 60000E-01  78  .130' 

8.0000  .12000  . 5370QE-0: 

100  100  . 20000E— 02  . 20000E-02 

PREHEATED  AIR  TEMP. 

FLAME  TEMPERATURE 

HEAT  EXCHANGER  OUTPUT  TEMP. 

HEAT  OUTPUT 

EXHAUST  OUTLET  TEMP. 

E.H.S.  EFFICIENCY 
COLD  START  PENALTY 
PRESSURE  DROP 
FLOW  FRICTION  POWER 
623.38  841.16  1388.4 


949. 9  K 

2374.9  K 
1096.3  K 
74.661  kw 

569.9  K 
87.675  V. 
87.361  gr 

2.193  KPa 
296.791  w 
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End  of  input  file 

*6 

6:*  100  100  0.002,0.002 
6s#  200  200  0.002,0.002 


B>EHS 

1.8040  10179.  V 

1.2500  . 40000E— 0 1  949.90 

.12000  . 60000E-01  61  . 

.18000  . 60000E— 0 1  78  . 

8.0000  .12000  .537001 

200  20q  . 20000E-02  . 20000E 

PREHEATED  AIR  TEMP. 

FLAME  TEMPERATURE 

HEAT  EXCHANGER  OUTPUT  TEMP. 

HEAT  OUTPUT 

EXHAUST  OUTLET  TEMP. 

E.H.S.  EFFICIENCY 
COLD  START  PENALTY 
PRESSURE  DROP 
FLOW  FRICTION  POWER 
624 . 00  837 .51  1 383 . 4 


&  3  AtHf 

10179.  V 


949.90  2.0000  1075.0  301.00 

61  . 13000E-02  3.8000  . 30000E-03  21.600 

78  . 13000E— 02  . 13000E-02  3.8000 
.  53700E-03 
20000E-02 

950.5  K 
2374.9  K 
TEMP.  1097.2  K 
74.613  kw 
570.4  K 
87.618  7. 

87. 187  gr 
2.179  KPa 
293.849  w 


3s  *  0.12,0.06,  61  1 . E-3, 3.8,0. 3E-3, 21.6, 

4:*  0.18,0.06,  78  1 . 3e-3, 1 . 3e-3, 3. 8, 


w. 

4:  * 

0.18, 0 . 06 , 

78  l.e— 

-s 

*E 

¥ 

B  >EHS  & 

% 

p 

1 . 8040 

10179.  V 

m" 

1 . 2500 

. 40000E— 01 

949.90 

h 

. 1 2000 

. 60000E-0 1 

61  U. 

. 18000 

.  60000E— 0 1 

78  .J. 

* « 

v. 

8 . 0000 

.  1 2000 

.  53700E 

V. 

• 

200  200 

. 20000E— 02 

.  20000E- 

2.0000 


1075.0 


301.00 


1 0000E-02  3.8000 


PREHEATED  AIR  TEMP. 

FLAME  TEMPERATURE 

HEAT  EXCHANGER  OUTPUT  TEMP. 

HEAT  OUTPUT 

EXHAUST  OUTLET  TEMP. 

E.H.S.  EFFICIENCY 
COLD  START  PENALTY 
PRESSURE  DROP 
FLOW  FRICTION  POWER 
645.59  818.03  1319.2 


976.9  K 
2382.0  K 
1082.8  K 
76.464  kw 
532.0  K 
89.792  •/. 
85.233  gr 
4. 493  KPa 
578.828  w 


ivs£*-; 


sicw 


.  •*.  i 

• *  j.  •  •  •  , 


s  V1 

/•M 


•W 
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List  of  Symbols  and  Abbreviations 

Mass  flow  rate 

Width 

Length 

Heat  flow 

Area 

Velocity 

Density 

Heat  of  vaporization 
Pressure 

Acceleration  of  gravity 

Shear  stress 

Radius 

Reynolds  number 
Drag  coefficient 

Angle  of  liquid  flow  to  horizontal 

Porosity 

Permeability 

Friction  coefficient 

Dynamic  drag  coefficient 

Position 

Profile  coefficient  for  momentum  flow 


Liquid 

Vapor 

Wick 


Hydraulic 


The  analysis  is  performed  on  a  finned  heat  pipe  evaporator  with  fins, 
length  width  b,  and  thickness  t,  as  shown  in  Figure  C-l. 

From  conservation  of  mass: 


£3 

L-  J 


From  conservation  of  energy: 


Combining  Equations  1  and  2  and  integrating  yields: 


(1  _  JL  ) 

K  a  '  x 


From  conservation  of  mass: 


m  =  p  A 


Combining  Equations  2  and  4  yields  an  expression  for  vapor  velocity: 


"  p  A  A 
v  v 


Combining  Equations  3  and  4  yields  an  expression  for  liquid  velocity: 
v.  .  <'-r>9 


p  z\x 
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,t in vyv^ w wjwv  .-v  a  A '.v.v.v.st.v.v.v. 


From  conservation  of  momentum,  the  steady  state  liquid  pressure  gradient  in 
the  evaporator  wick  in  the  direction  of  liquid  flow  is: 


dP„ 


2t, 


dx 


rh,  £ 


±  P  £g  sin  0 


(7) 


where  the  first  term  is  due  to  flow  resistance  and  the  second  term  is  due  to  gravity. 
Introducing  Reynolds  number  and  drag  coefficient; 


Re„ 


2rh,£  p£  vi 


2  tn 


pnY 


Combining  Equations  7,  8,  and  9  yields: 


dP„ 


dx 


-  -  (f„Re„  ) 


P£V£ 


'  '2rh  jl2  -  p£gSin 


Introduce  wick  porosity,  e  and  wick  area, 


A.  =  A  c 
£  w 


Define  wick  permeability,  K; 


_  2  enrh,e.a 

= 


Combining  Equations  6,  10,  and  12  yields: 


(8) 


(9) 


(10) 


(11) 


(12) 


S  -  .  '  .1 

•  ‘A 


•  *  «w*  V1 

m 

V  v\ 

KT*  *  V, 


V*  %v 


v;.*  i 


[.  1 


.  ■! 


V\  ‘ 

V  \! 


•  ..  •  *  vJ 


r  *m 

v>.  -n 


v* 


^£ 

dx 


K  AwX  p, 


P^g  sin  0 


(13) 
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Define  friction  coefficient,  F^;  and  substituting  into  Equation  13; 


F£  ~  KA  Ap. 
w 


(14) 


dx 


"F£Q  (  1  -  "X  )  '  P£g  sin  0 


(15) 


Integrating  Equation  15  yields  an  expression  for  the  liquid  pressure  drop 
(including  the  hydrostatic  term)  as  required  for  evaluating  the  capillary  pressure  in 
Equation  18  in  the  text. 

From  conservation  of  momentum,  the  vapor  pressure  drop  can  be  expressed 
as  follows: 


dP  -(f  Re  )  p  m  2m  dm 

V  V  V  V  V  0  V  V 

dy  ”  2A  r  7  p  ”  ^  A  2  p  dy 

’  V  V  V  V  V 


(16) 


The  construction  of  the  fin  is  such  that  the  vapor  path  is  packed  with  a  porous 
screen  structure  to  provide  structural  support  to  the  fin  walls.  This  support  structure 
is  analagous  to  the  screen  wick  used  for  the  liquid  flow,  with  regard  to  the  flow 
resistance  since  it  can  be  shown  that  the  vapor  flow  is  laminar  and  incompressible 
(see  Appendix  D).  Therefore,  by  introducing  the  vapor  screen  permeability,  Ky,  and 
vapor  screen  porosity  ey,  Equation  16  becomes; 


dP 
_ \ 

dy 
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Introducing  the  vapor  friction  coefficient,  Fy  and  the  vapor  drag  coefficient,  Dy  yields 


dP 
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dy 
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Integrating  Equation  18  yields  an  expression  for  the  vapor  pressure  drop  as  required 
for  evaluating  the  capillary  pressure  in  Equation  18  in  the  text. 
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